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SHAPE OF INITIAL PORTION OF BOUNDARY 
OF SUPERSONIC AXISYMMETRIC FREE JETS 
AT LARGE JET PRSSSURE RATIOS 
By Eugene S. Love and Louise P. Lee 


SUMMARY 


Calctilations have been made of the initial portion of the boundary 
of axisymmetric free jets exhaxisting at large ratios of jet static pres- 
sure to stream static pres stare fi^Dm a conically divergent nozzle having 
a jet exit Mach ntomher of 2-5 and a semidivergence angle of 15° • 
results of the calculations indicate the size and shape of the jet to 
be expected at large presstore ratios, the effects of the ratio of spe- 
cific heats, and the large Initial inclinations of the boundary that 
are likely to be encountered by hypersonic vehicles at hi^ altitude. 


INTRODUCTION 


In the proposed trajectories of most rocket -propelled hypersonic 
vehicles, the rocket propoilsion unit will remain in operation long 
enough for the exhausting jet to encoaonter the very low pressures of 
high altitudes. When this condition occurs, the ratio of the jet static 
pressure to stream static pressure becomes very large at the nozzle exit, 
and as a result the free jet expands greatly. The problems created by 
the presence of this lange, boolboris, free jet of gases immediately down- 
stream of the nozzle exit are, in the main, twofold. First, the large 
deflections of the free-stream flow caused by the jet flow may result 
in major aerodynamic interference upon nearby soorfaces; for example, 
large regions of separated flow may be created on the vehicle surface 
ahead of the rocket exit. Second, the heat existing in the core of 
this large mass of gases may at hypersonic speeds introduce direct heating 
problems through radiation by causing the temperatures of the nearby 
surfaces which are already experiencing high aerodynamic heating to 
exceed the critical values. None of the aerodynamic Interfeorence effects 
can be estimated with any reliability without some idea of the magnitoode 
of the region encompassed by the free jet and the shape of its bomdary. 
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For the purpose of gaining insight into the size and shape of free 
jets at large pressure ratios, a number of calculations have been made 
of the initial portion of the boundary of axisymmetric free jets 
exhausting into still air for jet static-pressure ratios ranging from 
about 60 to 42,000 for ratios of specific heat of 1 . 200 , 1.400, and 
1 . 667 . All calculations of the jet bovindaries were made for a conically 
divergent nozzle having a semidivergence angle of 15° and a Mach number 
of 2.5 for the jet flow at the nozzle surface immediately ahead of the 
exit. 


SYMBOLS 


The symbol notation that follows will be clarified by reference to 
the flow diagram shown in figure 1 . 


Mj jet Mach number; value at nozzle surface immediately ahead 

of exit 


Moo free -stream Mach n'umber; also Mach number ahead of exit 

shock or ahead of shock at beginning of separation 

p ambient or free-stream static presstire; also pressure ahead 

of exit shock or ahead of shock ~at beginning of separation 


jet static pressure; value at nozzle surface immediately 
ahead of eiit 


p' static pressure of external flow existing immediately ahead 

“ of jet exit and behind exit shock, or pressure in region 

of separated flow Immediately ahead of jet exit 

rj radius of jet exit 

X distance from plane of nozzle exit measiired parallel to 

jet axis 

y perpendicular distance from jet axis 

7 ratio of specific heat at constant pressure to -specific heat 

J at constant volume for jet flow 

5. initial inclination of free- jet boundary at exit 


Sjj semidivergence angle of conically divergent nozzle 
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CALCULATION OF BOUNDARIES 


General Procedure 

The method of characteristics was used throughout to calctilate the 
jet boimdaries. Leading characteristic lines (see fig. 1) and the initial 
values for the characteristic lines from the center of expansion were 
computed manually. Hie characteristic nets were determined on a card- 
programed electronic calculator. For expediency, only the boundary 
points and a portion of the characteristic net sufficient to define the 
Jet shock (by the foldback method of characteristics) were retained in 
the print -out phase of the automatic calculations. 

Although considerable experience had been obtained at the Langley 
Laboratory in previous calculations of Jet boundaries for low static- 
pressure ratios (of the order of 10 and lower ) , no information was avail- 
able for determining beforehand the increments between the characteristic 
lines Introduced from the center of expansion that should be used to 
Insure a reliable boundary solution for the large pressure ratios of this 
study. A similar situation existed as regards the spacing of the points 
along the leading characteristic line. After some trial calculations of 
boundaries, it was decided to introduce a new characteristic line from 
the center of expansion for approximately every degree of expansion. 

The spacing of the points along the leading characteristic line was 
selected on the basis of additional trial calculations and other factors, 
some of which are discussed subsequently in the section "Presentation 
of Boundaries . " However , it may be noted here that points along the 
leading characteristic line were calciilated for increments in 
of 0.02, and this spacing is recommended for the range of pressure ratios 
covered by this study. For values of and 0^^ greatly different 

from the values used herein ^Mj = 2. 

ratios greater than those covered in this study (of the order of 42,000), 
the rate at which the characteristic lines were introduced from the 
center of the expansion and the point spacing used along the leading 
characteristic line mi^t not be adequate. With regard to the latter, 
the shape and inclination of the leading characteristic line have, in 

general, first-order dependence upon and 0„. 

0 H 


= 15 ^) ^ and for pressure 


Static-Pressure Ratios for Calculations 

The choice of static-pressure ratios for the calculations was 
influenced to some extent by an examination of the variation of the ini- 
tial inclination of the Jet boundary 6j (see fig. 1) with static- 

pressure ratio Pj/Pjj^. This variation is presented in figure 2 for 
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the values of the ratio of specific heats used in the boundary cal- 

culations. It is interesting to note first the magnitude of 5j for 
the nozzle with a sonic exit = 1, “ ®°) f since at very large 

values of this nozzle gives the largest possible values of 6j 

(fig. 2) . By comparison with these values for the sonic nozzlej the 
values of Sj for the nozzle of this study ^Mj = 2 .^, “ ^5°) 

significantly lower at very large values of Pj^p^, and this difference 

increases with decreasing Of course at very low values of Pj^Poo 

any conically divergent nozzle will produce a value of Sj greater than 

that for the sonic nozzle, since 5j = 0jj at Pj^P„ = 1; this 
instance, 5j = 0jj = 15°- Pj/P« = !■ 

The curves of figure 2 Illustrate that, at large values of Pj^p^, 

conically divergent nozzles likely to be encountered in practical applica- 
tions may produce large initial inclinations of the boundary. For rockets 
where y^ is of the order of 1.2 and sometimes less, this initial incli- 
nation may easily exceed 90°, and at extremely large pressure ratios it 
is not inconceivable for the free- jet flow at the jet boundary to be 
directed initially alinast 180° opposite to the direction of flow along 
the jet axis. Thus, to the general objective of gaining some familiarity 
with the size and shape of free jets at large pressure ratios is a^ded 
the objective of calculating some boundaries for which 6j > 90°. In 

addition, boundaries for which a comparison will show the effects of y ^ 
for similar values of and for similar values of Bj are desired. 

The curves of figure 2 show that the variation of 6j with 

is such that at Pj^Poo order of 10® for y^ = I.667 and of the 

order of ID^^ for 7^ = 1.200, the initial inclination has approached 

closely its maximum value. Consequently, if calculations could be made 
near these values of Pj/Poo (°^ higher), they shoTild be sufficient for 

most purposes to represent jets operating in the range from these values 
to infinity. Unfort\mately, altho-u^ solutions were attempted for these 
extremely large press\rr'e ratios, satisfactory resailts could not be 
obtained without successive refinements and reprograming of the auto- 
matic calculations. For these reasons and because of other commitments 
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of the electronic calculators, the boundaries at these extremely large 
pressure ratios vere not obtained. Bie values of and the cor- 

responding values of 6j, for which boundaries were calculated are given 
in the following table ^boundaries for order of 10 and 

less were available from previous calciolatlons) : 


= 1.667 

7 = 1.400 

7 = 1.200 

j 


deg 


6j, deg 


5 j, deg 

59.92 

569.2 

l,.6l8 

lt 8.55 

57.72 

6o.7li- 

81.94 

1,346 

2,692 

60.02 

75.59 

78.55 

76.58 

1,303 

9,776 

41,820 

70.99 

94.51 

107.6 

115.6 


The exact values of Pj/Poo ^sed in the calculations were determined by 

the arbitrary condition of terminating the expansion at one of the 
manually calculated characteristic lines from the center of expansion. 


RESULTS AND DISCUSSION 


Presentation of Boundaries 


The calculated boundaii.es are presented in the form of curves with 
discontinuous slopes } this type of presentation is compatible with the 
lattice-point characteristic method (that is, calculated points on 
boundary joined by strai^t lines). In practical applications a smooth 
curve is, of coxirse, fitted to the die continuous -slope form of the bound- 
particular care being taken to maintain the Initial inclination of 
the boundary. 


Separate presentation with jet shocks .- The calculated boundaries 
together with the indicated position of the shock formed within the jet 
(given by the intersection and coalescence of characteristic lines of 
the same family) are shown separately for each value of -p in fig- 


Is 


ures 3 to 5 * 7 j = 1 . 200 , only the boundary for pyp^= 76.58 

presented separately (fig. 5) inasmuch as only a minute segment of the 
jet shock was Indicated in the inadvertently shortened calculation at 
p ^p = 1,503, and none was indicated at the higher values of Pj^P^* 
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All bo\mdarles for 


7 j = 1.200, as well as for the other values of 




are presented together subsequently. At first glance the results for 
all values of 7 j give indication of a tendency for the jet shock to 

be eliminated with increasing Pj^P large Pj^P • however, 

appear debatable as to whether a valid indication of such an occurrence 
could be substantiated by these calculations, inasmuch as an inherent 
stretching of the meshes within the characteristic net occurs with 
increasing Pj^P^o these calc\ilations . Nevei-theless, it may be rea- 
soned that for Pj^Poo “ “ velocity at the boundary is radial; that 

is, the characteristic line is in the same direction as the direction 
of the velocity, which is the direction of the boundary. Consequently, 
the only coalescence of the characteristic, lines occurs on the boundary. 

It follows that if at Pj/Poo “ “ this coalescence were indicative of a 

shock, the shock would coincide with the boundary. However, since a 
shock nrust have a pressure rise across it, the existence of a shock on 
the boundary represents an impossible condition for Pj/Poo “ “ ^ view 

of the following; first, a pressure rise cannot exist across a boundary, 
and second, the existence of a pressure rise at Pj^Pco ~ “ would lit^jly 

that, in expanding to a vacuxnn, the Jet flow has expanded to a pressure 
less than a vacuum and thus requires a pressure rise at the boundary to 
restore its pressure to the ambient vacuum, a clearly impossible situa- 
tion. Thtis, it appears that at Pj^P^ = “ no jet shock will be present, 

and an indication of the elimination of the jet shock with increasing 
Pj^p^ may, at large 




be reasonable. 


In figures 4(a) and 5^ "two boundaries and two shocks are shown. 

In both figures these two boundaries and shocks represent the examina- 
tion of the effect of point spacing along the leading characteristic 
line upon the shape and position of the calculated boundary and the jet 
shock. QMs examination was primarily an outgrowth of the desire to 
retain a programing capacity (50 points along the leading characteristic) 
for the automatic calculations that had been used in previous stiidles 
at lower pressure ratios and, at the same time, to obtain, if possible, 
a significant portion of the boundary and jet shock without inciirring 
large errors. Ihe shorter boundaries and jet shocks correspond to cal- 
culations made by the -use of the first 50 points along the leading c]mr- 
acterlstlc ^that is. Increments in maintained standard at 0.02^. 

Ihe longer boundaries and jet shocks correspond to results obtained by 
use of every other point of the first 50 calculated points on the leading 
characteristic ^increments in of 0.04J, every third point of the 

second 50 ^increments in x^Tj of 0 . 06 ^, and every fourth point of an 
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additional 20 ^increments in 0.08^. The results show that 

the errors associated with the larger point spacing are small; conse- 
quently this spacing was used in most of the calculations. In these 
figures and those to follow, a few of the boundaries do not have the 
extent to be expected from this point spacing. For these boundaries, 
difficulties encountered in the automatic phase of the calculation 
shortened the calculations. Since with one or two exceptions the com- 
pleted portions of the boundaries for these shortened calculations are 
believed to convey adequately the shape of the boundary over a consider- 
able distance from the jet exit, no effort was made to continue these 
calculations . 

Composite presentation .- In figure 6, the boundaries for each value 
of 7 j are presented compos itely to show more readily the effect of 

increasing Pj^p^^. Included are boundaries for Pj^Poo “ ^ cal- 

culated in a previous study. The enormoTis size that the free jet may 
attain at large pressure ratios is readily apparent, particularly at 
7 j = 1.200 (fig. 6(c)) when the initial inclination of the boundary 

approaches and exceeds the vertical. 


Effect of 7 j 


Effect of 7 j for essentially equal pressure ratios.- Some idea 
of the effect of 7 ^ at large pressture ratios upon the shape of the 
boundary for essentially equal values of Pjy^oo obtained by 

comparing the bovLndarles for 7 j = 1.400 and 1.200 at Pj/Peo *** ^>500 as 
shown in figure 7- (Actually, for 7 ^ = 1.400, Pj/Poo ~ 1^546 and for 
7j = 1.200, Pj/p„ = 1 , 305 .) In spite of the short extent of the 
boundary for 7 ^ = 1 . 200 , the comparison shown in figure 7 is sufficient 
to point up the large effect that has upon bOTindary shape at large 

pressure ratios. 

Effect of 7 j for essentially equal Initial inclination.- At low 
values of Pj^Pco order of 10) the procedure of simulating a 

jet of one value of 7 j by a jet of another value of 7 j thro'ugh dupli- 
cation of the jet inclination of 5j by changing the value of Pj^^Poo 

is a useful artifice in many experimental investigations. An examination 
of the adequacy of this proced-ure at large pressure ratios is of interest. 
Figure 8 compares the boundary of the jet for y^ = 1.400 with the 
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boundary for = I.667, both having an initial Inclination of about 6o° 

but widely different values of Pj^p^- (The value of 

7j = 1.667 is about 20 times that for y^ = 1.400.) The comparison is 

sufficient to show that this procedure does not yield accurate simulation 
of the boundary at large pressure ratios and thus cannot be expected to 
give other than first-order simulation of interference effects . For 
exrbremely large values of Pjy^Poo curves of figure 2 serve to demon- 

strate that simulation by this procedure is impossible for the case of 
primary interest, namely, simulation of a jet having low 7j by a jet 

having a higher 7j (with significant differences in y^, of course, say 

of the order of 0.1 or more) . For example, it is inpossible for the jet 
of this study with y^ = 1.400 to duplicate the value of for the 

jet with 7j = 1.200 when the latter jet operates at values of Pj^p^ 

greater than about 10^. 


Comments on Bo-undaries at Large Pressure Ratios 
With Supersonic Ambient Stream 

I&pothetical hypersonic vehicle ,- Althou^ the calculated boundaries 
of this study are for the jet exhatisting into still air, it is apparent 
from these calculations that a rocket-propelled hypersonic vehicle may 
encounter serious jet interference problems if it encounters the high 
pressure ratios that are likely to occur at high altitudes. Some comments 
on the jet issuing from a hypothetical ground-launched hypersonic vehicle 
are thus considered in order. 

Figure 9 presents an ass-umed speed-altlttide variation of a hypo- 
thetical hypersonic vehicle. For simplicity a single rocket is assumed 
to provide power for the vehicle and the rocket is assumed to operate 
at design value at sea level ^that is, at Pj/Poo = ^ practical 

applications the value of Pj^oo usually less than 1.0 

if the rocket is to encounter high altitudes. 

Variation of pressure ratio with altitude .- The variation of Pj^Poo 

with altitude was calculated with p assumed equal to the pressure at 

00 

altitude and with the pressure at altitude varying according to the 
ICAO standard atmosphere (ref. 1) with extensions to higher altitudes. 

The value of p. was assumed to be invariant with altitude. The 
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variation thus obtained is shown by the solid curve in figure 10 and 
represents the conditions shown in the following sketch; 



It is iir^jortant to note that the solid-line curve in figure 10 is inde- 
pendent of the value of at altitude; the speed of the ambient flow 

could be zero at all altitudes or covild vary in some arbitrary manner 
as in figure 9, 

In order to examine the effect of speed-altitude variation of the 
hypothetical vehicle (fig. 9) on the initial inclination of the ;}et 
boundary consideration must be given to the pressure rise throxigh 

the shock in the ambient flow. For large pressure ratios separation 
will undoubtedly occur on the external surface ahead of the base 
because the boundary layer will not withstand a pressure rise of the 
magnitude necessary to turn the flow throu^ the large initial deflec- 
tions of the Jet boundary. The following sketch probably represents 
a more realistic representation of the phenomena; 
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Little is known at present about the pressure rise associated with 
separation at hypersonic speeds ; consequently two conditions have been 
examined. For one condition a pressure rise corresponding to 

a 50° turning of the flow through an oblique shock at separation is 
SLSSumed to occur; for the other condition, in order to have a quasi 
limiting condition, the press-ure rise through a normal shock is assumed 
to have taken place ahead of the separated region. In the determina- 
tion of the restO-tlng values of the pressiire p'^o is assumed 

to hold throughout the separated region. The curves corresponding 
to the oblique shock with 50° turning and to the normal shock are shown 
in figure 10 . Of course the portions of the curves at very low values 
of Pjy^p'^ have no meaning since a 50° deflection of the external flow 

by the jet would be impossible at these low values. For large pressure 
ratios, however, the oblique-shock curve shotild give a reasonable first- 
order indication of the values of pj^p'^ that would be encountered, 

and the fact that the normal-shock curve is not greatly removed from the 
oblique -shock cvirve is believed to indicate that values of 5 ^ calculated 

from the oblique-shock curve could not be greatly different from those 
that are experienced with practical configurations. 


Variation of initial bo-undary inclination with altitude.- Since 


the curves of figure 10 are independent of 7j, 6jj, and Mj. 

M. - 2.5, e, 


be applied to the nozzle of this study having 
varying values of 7^ . 


they may 

15°, and 


j TO 

Values of the initial inclination of the jet 


boundary 6j for this nozzle have been determined with the aid of fig- 
ures 2 and 10 and the res-ults are shown in figure 11 . The solid-line 
curves are for = 0, and the dashed-line curves are for the speed- 


altitude variation of figure 9 with the oblique-shock separation of the 
bo-undary layer. Although only the curves for 7j = 1.200 may be con- 
sidered applicable to the irocket-propulslon -unit of the hypothetical 
vehicle, the curves for the other values of 7. have been included -to 

U 

point up the effect of 7. and for possible information on the behavior 


of jets on reaction controls. Again the lower portions of the dashed- 
line curves ^ small 

solid-line curve and -the dashed-line curve at a particular 7j shows 

the effect that the presence of supersonic ambient flow, defined 
according to figure 9 -with separation ahead of the base, has in 
reducing the initial deflection of the boundary. As the extremely large 
pressure ratios of high altlt-ude are approached (fig. 11), the effects 
of the ambient flow may be readily deduced from -these c-urves to be 


6j^ have no meaning. The difference b e-tween the 
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secondary In determining the value of 



Of particular importance is 


the indication that feasible rocket-propelled hypersonic vehicles may 
encotinter within their range of operation Initial inclinations of the 
boundary approaching and exceeding 90 °; consecjuently, large regions of 
separated flow, together with large aerodynamic and possible aerothermal 
Interference effects, may occur. 


Remarks on boundaiy shape .- iiie indication that the effects of the 
ambient flow are secondary at extremely large pressure ratios in deter- 
mining the value of 5j does not imply that the shape of the boTindary 

will not be significantly affected by the ambient flow at lower pressure 
ratios. No reason is apparent why the procedure employed at low pres- 
sure ratios for determining the shape of the portion of the jet boundary 
that bounds the region of separation (from A to B in sketch (b)) should 
not be useful at large pressure ratios; in this procedure this portion 
of the boundary is determined by employing the boundary of the jet 
exhausting into still air at a pressure ratio equal to Pj^p'„. However, 

at point B of sketch (b) , where the separated region ends and the jet 
flow and ambient flow interact directly, a large and abrupt turning 
of the boundary may be expected at large pressure ratios and hypersonic 
speeds, if the flow deflections accompanying separation at these speeds 
are of the order of 30 ° (turbulent boundary layer preceding s^aration), 
as seems likely on the basis of available experimental information. 

Beyond point B the ambient flow will have significant influence on the 
boundary shape in the same manner as it does at low pressure raid. os and 
supersonic speeds. 


CONCLUDING REMARKS 


Calculations have been made of the initial portion of the boundary 
of axlsymmetrlc free jets exhausting into still air when the ratio of 
jet static pressure to stream static pressure is high . Hhe jet Mach 
number of these caJLculatlons was 2.5 and the semidivergence angle of 
the nozzle was 15 . Die calculations covered static-pressure ratios 
from about 6 o to 42,000 and ratios of specific heats of 1 . 200 , 1.400, 
and 1 . 667 . T3ie results indicate the enormous size that the free jet 
may attain at large pressure ratios and point up the large effects of 
the ratio of specific heats at large pressure ratios. A cxursory examina- 
tion of the behavior of the propulsive jet on a hypothetical rocket - 
propelled hypersonic vehicle Indicates that feas-ible hypersonic vehicles 
may encounter within their range of operation initial Inclinations of 
the free- jet boxmdary approaching and exceeding 90° • Consequently, 
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large regions of separated flow may occvir on such vehicles, together 
with large aerodynamic and possible aerothermal interference effects . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 14, 1957* 
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Figure 1.- Slcetch of jet exhausting from conically divergent nozzle into 
atm air shyowlng pertinent details nrui parameterB. 
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(a) -i= 81.94; = 60.02°. 


figure 4.- Calculated boundaries and jet shocks for jet exhausting into still air with 

7 j * 1.400, Mj => 2-5, and 6^ » 15°. Longer boundary and shock correspond to greater spacing 

between points on leading characteristic theai for shorter boundary and shock (see text). 


56m VDW 



NA.CA TN 1<-195 


19 



Figure I)-.- Continued. 
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(b) 7j = 1.400. 


Flgiore 6.- Continued. 
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7j = 1.200. 
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Figure 9 .- Assumed Mach number altitude variation o 

vehicle . 




if hypothetical 
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= 1 at sea level, pj is invariant with altitude, and p„ 


equals pressure at altlt\ide. 
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Flgiire 11.- Variation, of initial Inclination of boundary with altitude. 

Mj = 2.5J 9if = 15°. 


NACA - Langley Field, Vj.. 


